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Abstract 

The influence of surface cooling on compressible boundary-layer instability is dis- 
cussed theoretically for both viscous and inviscid modes, at high Reynolds numbers. 
The cooling enhances the surface heat transfer and shear stress, creating a high- 
heat-transfer sublayer. This has the effect of distorting and accentuating the viscous 
Tollmien-Schlichting modes to such an extent that their spatial growth rates become 
comparable with, and can even exceed, the growth rates of inviscid modes, including 
those found previously. This is for moderate cooling, and it applies at any Mach num- 
ber. In addition, the moderate cooling destabilizes otherwise stable viscous or inviscid 
modes, in particular triggering outward-traveling waves at the edge of the boundary 
layer in the supersonic regime. Severe cooling is also discussed as it brings compressible 
dynamics directly into play within the viscous sublayer. All the new cooled modes found 
involve the heat- transfer sublayer quite actively, and they are often multi-structured in 
form and may be distinct from those observed in previous computational and experi- 
mental investigations. The corresponding nonlinear processes are also pointed out with 
regard to transition in the cooled compressible boundary layer. Finally, comparisons 
with Lysenko and Maslov’s (1984) experiments on surface cooling are presented. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 


Surface cooling is fairly often applied on supersonic-hypersonic flight vehicles, and in 
some cases it is essential for the preservation of the vehicle itself, especially at higher Mach 
numbers. The effects of such cooling on the instability and transition properties of the local 
compressible boundary-layer flow are then of much concern. Given the practical importance, 
it is perhaps surprising that there are not more experimental, theoretical and computational 
studies of surface-cooling effects on instability and transition; see the experiments of Van 
Driest and McCauley (1957), Wisniewski and Jack (1960), Richards and Stollery (1966), 
Stetson and Rushton (1967), Van Driest and Blumer (1968), Kendall (1975), Demetriades 
(1978), Lysenko and Maslov (1981,1984), and the theoretical suggestions and computations 
of Lees and Lin (see Lin 1955), Reshotko (1963), Mack (1969,1975,1984,1986), Maslov (1974), 
Gapanov and Maslov (1980), Malik (1987). The flow phenomena produced can be complex, 
and the findings and suggestions from the works above are a little mixed, but Lysenko and 
Maslov (1984) seem to provide a reasonable overall summary, namely that the so-called first 
modes tend to be stabilized by surface cooling whereas the second modes are destabilized to 
a limited extent. Also, islands of first-mode instability have been predicted by parallel-flow 
theory at the lower Reynolds numbers, although these suggestions are questionable due to 
the nonparallel-flow effects present then: see also later. 

The present theoretical investigation, in comparison, suggests that the effects of surface 
cooling can be much more pronounced and “dangerous.” In particular it is found that the 
cooling, even if only moderate, so distorts the viscous-inviscid Tollmien-Schlichting (TS) 
modes [first modes] that their spatial growth rates can become comparable with or larger 
than those of the inviscid modes [second and higher modes] described above, at high Reynolds 
numbers. This reversal of the usual ordering of the instabilities arises because of the enhanced 
heat transfers and shear stresses produced in the basic boundary-layer flow at the surface 
by the presence of cooling, which results in a sublayer of high heat transfer emerging on the 
surface. The increased surface shear stresses cause the typical viscous-inviscid TS wavelength 
to be decreased and the flow structure associated with the TS wave is soon altered, as the 
surface is cooled, bringing about the rise in growth rates, as well as rendering previously 
stable modes unstable. 

The present predictions appear to agree broadly with Lysenko and Maslov’s (1984) sum- 
mary finding above, and with their experiments for supersonic boundary layers (see also 
the comparisons in §6 below), in that the wavenumbers and frequencies associated with the 
boundary-layer’s first-mode instability increase with surface cooling, at a given Reynolds 
number. There is however some possible disagreement as well (which may be a matter of 
interpretation) in the sense that the theoretical growth rates increase also, as described in 
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the previous paragraph. Thus the prediction here is that, with cooling, instability may be 
delayed but when it occurs it is likely to be more violent than in the non-cooled, e.g. insu- 
lated, case. This is mainly for the viscous-inviscid first modes, but similar findings emerge 
for the inviscid second modes also as discussed later. 

“Moderate cooling” of the surface, and its influence on the compressible viscous modes, is 
dealt with in §§2 and 3 below, based on the theory in Smith (1987). The new flow structure 
produced is examined first, followed by its instability properties. The flow structure is of 
the compressible Rayleigh inviscid kind across the majority of the boundary layer, but in a 
quasi-steady form, which admits pressure-displacement interaction with the viscous sublayer 
response. Unlike the non-cooled TS instability modes, which in supersonic flow must be 
directed outside the wave-Mach cone, i.e. be effectively subsonic, the present surface-cooled 
ones can be unstable for any obliqueness angle, including the two-dimensional (2D) case. 
This and other factors of the link between the cooled and the non-cooled TS modes are also 
explored in § 3 , along with the influence of cooling the surface beyond the moderate stage. 
There the spatial growth rates become still more accentuated, for any Mach number, as the 
viscous-inviscid disturbance is concentrated relatively near the surface, and flow-structural 
change re-occurs. We devote a separate section, §4, to the resulting stage of “severe cooling” , 
to emphasize its novel nature. This concerns principally the intrusion of compressibility 
effects into the dynamics of the thin viscous sublayer. Although the work overall is mostly 
on viscous modes, pure inviscid modes are addressed in §5 as a result of the behavior found 
earlier in the moderate-cooling regime. Once again, cooled-surface modes which appear to 
be distinct from those found in earlier studies are implied, for any Mach number with 
the heat-transfer sublayer playing an active part throughout. The new inviscid modes, like 
the viscous ones of §§ 2 - 4 , all have relatively small phase speeds, and they include some which 
radiate, i.e. provoke outward-traveling waves at the edge of the boundary layer. 

Nondimensional variables are used such that the freestream velocity is ( u , v, w ) = (1,0, 0), 
in cartesian coordinates (x, y, z) where the airfoil chord is unity. The nondimensional pres- 
sure, density, temperature and viscosity are written p,p,T,p. in turn, equal to p,*,, 1,1,1 in 
the freestream, where Poo = 1/(7 Af^), and t denotes the usual convectively scaled time. 
The basic boundary layer then has 0(Re~ lf2 ) characteristic thickness where Re is the global 
Reynolds number produced above and the boundary-layer properties are defined by the pro- 
files ( u,T,p ) = (U 0 ,T 0 ,po)(y), with y = Re~ 1/2 y, at the local station x = x 0 , z = z 0 under 
consideration. In §§ 2-4 nonlinear viscous-inviscid interactions are considered first since the 
understanding of flow transition is the ultimate theoretical goal, although linear solution 
properties are then treated, as a first study, since this type of approach for surface cooling 
appears to be a fresh one. The corresponding nonlinear inviscid versions in §5 can likev.ise 
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be set up. Further discussion is presented in §6, along with some comparisons with the 
experiments of Lysenko and Maslov (1984) and related comments. 


2. THE MODERATE-COOLING REGIME 

For order-one values of the temperature ratio T^/T^(= T w ) the compressible three- 
dimensional TS modes are controlled by the triple-deck structure as described by Smith 
(1987), as long as the other parameters involved are not too extreme. That structure pro- 
vides the springboard for the present theory on surface cooling and so its underlying form is 
noted here. The expressions for the flow solution in the three decks are: 


upper deck, y = Re 3/8 Kim 1/2 y u and 
[u,v,w,p] = [1,0,0, + i7e _1/4 C 1/4 A 1/2 m _1 / 4 [t^ 2 ),m 1 / 2 ^ 2 ), : u^ 2 ),p( 2 )]; J 
main deck, y = Re~ x ^ 2 y and 

[u,v y w t p] = {U 0 ,b,0,p oo } + 0[Re- 1 l*,Re- 1 l A ,Re- x l A ,Re- x l i ]\ J 
lower deck, y — Re~ s ^ 8 m~ 1 ^ 8 C 5 ^ 8 T 3 ^ 2 X~ 3 ^ A Y* and 

[u, v, w,p] — [Re~ 1/8 C 1/8 T 3/2 X 1/4 m~ 1/8 U*,Re~ 3/8 C 3/8 T^/ 2 X 3/A m 1 ^ 8 V*, 


(2.1a) 


(2.16) 


(2.1c) 


i2e - 1 /8 c 1 / 8 jU/2 A l/ 4 m - 1 /8^* i2e - 1 / 4 C7 1 / 4 A 1 /2 m - 1 / 4 p,]; 

throughout, the streamwise and spanwise length scales and the time scale are given by 


(x-x 0 ,z- z 0 ) = Re-^K^X^Z*), t = Re-^ 4 m-^ 4 T w C l/4 X~ 3/2 t* . ( 2.1 d) 

Here m = | — 1|, and P and — A are the unknown pressure and displacement variations 
(independent of y,Y*),C is the Chapman constant, A(x 0 , *o) is the 0(1) skin-friction coef- 
ficient of the incident compressible boundary layer and (xo, zo) is the local position under 
consideration, while K\ = C 3 ^ 8 T^/ 2 A -8 / 4 m 3 / 8 . Similar expansions hold for the temperature 
T, the density p and the viscosity p. As a result the overall nonlinear three-dimensional TS 
behavior here is governed by the lower-deck equations and boundary conditions, 


U* x . + Vf. + W* z . = 0 , (2.2a) 

u;. + u*u* x , + v*u Y . + w*u* z , = -p* x .(x\ z*, t*) + u^. Y ,, ( 2.2 b) 

W* + U*W* X . + V'Wl, + W*W' Z . = -P' z .{X\Z\t*) + (2.2c) 

u* = V* = W* = 0 at Y* = 0, (2. 2d) 

IT + A\X\Z\f), W* -+ 0, as Y* -> oo, (2.2e) 
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coupled through pressure-displacement interaction with the upper-deck linearized potential 
flow problem (for supersonic and subsonic flow respectively) 

[(Mi - l)(d\, T dl) - aj.]p« 2 > = 0, (2.2/) 

where 

-> 0 as y u -» oo; -► P*, pj, 2) -► A* x * x . as y u -► 0 + . (2.2 g) 

Suitable constraints at large \X*\ } \Z*\ } are also assumed, and imposed forcing can be incor- 
porated in either or both of the parts (2,2a-e), (2.2f,g). Various nonlinear features of the 
unsteady viscous-inviscid system (2.2) for nonlinear TS evolution are addressed by Smith Sz 
Walton (1989), Smith (1988, 1989) among others, while the linearized version corresponding 
to relatively small compressible TS waves is studied by Smith (1987). Here we follow the 
same procedure of setting up the nonlinear stage first and linearizing afterwards, in dealing 
with surface-cooling effects. 

The main impact of surface cooling, i.e. small T w , is felt partly through the T w factors 
themselves in the scales above and partly through the alterations in the mean-flow profiles 
Uoy Tq near the surface, specifically in the large heat transfer and skin friction produced there. 
This can be seen in two ways. One is by means of the Howarth-Dorodnitsyn transformation, 
y' = J o r 0 - l dfc since U 0 ,T 0 are regular in terms of the transformed variable y'. So if T w is 
reduced then initially T 0 oc y' at small y' with an 0(1) constant of proportionality. Hence 
the transformation implies that y oc f 0 T 0 dy r is quadratic in y\ and therefore the temperature 
profile Tq is tending to zero like the square-root of the normal distance y. The same applies to 
the velocity f/ 0 . That implies large temperature and velocity gradients developing near the 
surface. The precise size of these is determined by a sublayer produced where y is O(T^), or 
y 1 ~ and the low surface temperature becomes comparable with the linearly decreasing 
(in terms of y f ) contribution above, so that T 0 « T w + y' times 0(1), i.e. T w + 0(1) times 
[y + O(l)] 1 / 2 . It follows that the surface heat transfer, dT 0 /dy at y = 0, is 0{T~ l ). Likewise, 
the surface shear stress is 0{T~ l ). The alternative is to consider the basic-flow momentum 
and heat-flux equations, which are dominated by their viscous contributions near the surface 
at small y, i.e. by d(fiodUo/dy)/ dy = 0 and 8{fi 0 dTo/dy)/dy = 0 respectively. Since /x 0 oc T 0 
for a Chapman fluid (see also later) the second balance above integrates to give T 0 dT 0 /dy 
uniform in y and hence Tq is linear in y; so again the almost square-root response in the 
temperature T 0 and the large 0(Tj J ) heat transfer 8T 0 /dy are implied at small y and T w . 
Similar results for Uq stem from the momentum balance above. The effect of surface cooling, 
then, is to induce a new sublayer in which, for certain 0(1) constants , c 2 , 

y = T 2 J, To = T w {\ + 2 dy) 1 / 2 , Uq = T w —{(1 + 2 Cl yf' 2 - 1}, (2 .‘6a) 
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and in consequence the scaled heat transfer and skin friction 


(dTo/dy) 5=0 ~ T~'c u X = ( dU 0 /dy )y= 0 ~ T~ l c 2 , (2.36) 


are increased substantially. 

The implications for the TS waves follow directly from the increase in A in (2.3b) and the 
decrease in T w . Thus the TS length scales of ( 2 . 1 a,c,d) all shrink now, as might be expected 
physically, and in particular the streamwise length and the maximum normal distance both 
decrease like T^J 2 . Similar deformations occur in the typical lower-deck thickness and in 
all the velocity and pressure scales, for instance, and eventually a new flow structure must 
come into play as T w continues to drop. Various options present themselves but the first 
distinct new structure is found to arise when T w is decreased to the order iie -1 / 12 , at which 
stage the upper-deck extent becomes coincident with the main boundary layer’s, as does 
the streamwise length scale. So novel elliptic effects must then enter the boundary-layer 
dynamics. 

In the new stage of “moderate cooling”, therefore, we have 

T w = Re-^ 12 f w , (2.4) 

say, with T w typically of order unity, and a three-tier structure comes into operation. See 
Figure 1. The scales and balances involved are suggested by the ones mentioned earlier but 
combined with (2.4). Across most of the boundary layer, the main tier where y is 0(1) is 
given by 


[u, V, w, p, p, T ] = [U 0 , 0, 0, Poo, Po, To] + Re 1/4 pX 2 c 2 3 [u, v, w, p, p, T] + ■ ■ ■ (2.5) 

with c 3 = {C /Xp 2 ) 1 / 3 , while in the buffer tier where y = Re~ 2 ! 3 y 

[u,v,w,p, p,T] = [f2e _1 ^ 12 {t + Re^^c^Au' , — Re^^c^Axu, Re~ l ^ 4 X 2 c\D /u, 

Poo + Re~ l/4 pX 2 clP, Re 1/12 p, Re~ 1/12 T W ], (2.6) 

to leading order, with Dx — ~Pz, and in the lower tier y — Re~ 3 ^ 4 c 3 Y with 

[u,v,w,p,p,T] = [Re- 1/6 Xc 3 U,Re-* /12 Xc 2 3 V,Re- 1/6 Xc 3 W, 

Poo + Re-^ 4 pX 2 c\P, Re^p, Re-^ U T W }. (2.7) 


Here the extra factors C, A(= c 2 /T w ), p(= T^ 1 ) present are all 0(1) and are included to pro- 
duce a normalized equation set below, taking account of the scaled surface density, viscosity 
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and shear values which would otherwise appear in the lower-tier system for example. The 
streamwise and spanwise length and the time scales now are defined by 

(x — x 0 ,z — z Q ) = Re~ 1 ^ 2 (X, Z), i=— = . (2.8) 

AC3 

Hence the main-tier properties are controlled by the linearized quasi-steady inviscid balances 


Po u x + PxUq + PoVy + po v v + powz = 0, (2.9a) 

U Q u x + vU 0v = - Px/po , (2.96) 

Uo(v x , w x ) = -( py , Pz)/p 0 , (2.9c) 

Uo(poPx - 'yPooPx) = IPooPayV, (2.9 d) 

leading to the compressible pressure equation 

Pyy + Pzz + C 1 - M o )Pxx = 2 Mo-p-/ M 0) (2.10a) 


with effectively zero phase speed. Here M 0 = UoM^/Tq^ 2 . The constraints on p are 


p —* 0 [or outgoing waves] as y — > oo, (2.106) 

2 

P ~ P + g rA XX y 3/ 2 as j/ -> 0+, (2.10c) 

together with appropriate boundedness conditions far upstream and downstream (see below), 
where (2.10b) reflects the match with the oncoming undisturbed free-stream motion and 
(2.10c) is necessary for merging with the buffer tier. The 3/2-dependence in (2.10c) comes 
from analysis of (2.10a), given that the basic-flow profiles Uo,Tq exhibit the square-root 
dependence 

C/ 0 ~c 2 (2/c 1 ) 1/ Y /2 ) T q ~ (2c!) 1/2 y 1/2 , as y -> 0+, (2-11) 


consistent with (2.3). The 0(1) parameter r introduced in (2.10c) is defined by 

= <j /3 2 l ' 2 Tl 
7 c\ l2 C l l z ’ 


(2.12a) 


for consistency with the buffer tier below. Again, the associated velocity perturbations 
respond in the form 

(u, v, w) ~ (y _1/2 , y 1/2 , 0(1)), (2.126) 


at small y , agreeing both with (2.10c) and (2.6). The buffer tier of (2.6) is also governed by 
predominantly inviscid linear dynamics but it accommodates the basic inner-layer behavior 
(2.3a), thereby smoothing out the 3/2-dependence, since the velocity profile 

6(5) = -[(?’ + 2 c.5) 1 '* - f.J, (2.13) 

C 1 
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is regular. The final, lower, viscous tier of (2.7) then yields the three-dimensional unsteady 
nonlinear bound ary- layer equations, 


jj x + Vy + Wz — 0, (2.14a) 

Ui + UUx + VU Y + WUz = - Px{X , Z , i ) + Uyy , ( 2 . 146 ) 

Wi + UWx + VWy + WWz = — *) + Wyy, ( 2 . 14 c ) 

subject to the conditions 

U = V = W = 0 &tY = 0, (2.1 4d) 

t7~y + A(X,Z,i), W — ► 0, as Y — ► oo, (2.14e) 


for no slip at the surface and to match with the buffer tier’s response in turn. 

The nonlinear system describing moderate cooling is therefore given by the viscous- 
inviscid interaction of (2.10a-c) combined with (2.14a-e). A match may be verified with the 
earlier structure of (2.1) at large t, corresponding to large T w , and certain extra aspects of 
this are examined in the next section. Of most concern here are the properties of the new 
cooled-surface TS system at 0(1) values of r and the effects of the inviscid but rotational 
elliptic part (2.10) of the interaction, which represents the major new feature present in the 
moderate-cooling regime. Further, small-disturbance properties can be of much interest, for 
practical and other reasons including the fact that the length scales and hence the spatial 
growth rates expected are now comparable with those of pure inviscid modes, and so we 
address here the linearized instability version in which (77, V, W, P, A , p) are taken to be small 
perturbations of the basic undisturbed state (Y, 0,0, 0,0,0). This has the effect of keeping 
(2.10a-c) intact but with the lower-tier response becoming linearized. In consequence cool- 
surface TS waves oc exp(i(aX + j3Z — Q, i)) are possible with wavenumbers a, (3 and frequency 
ft, in which case we are left with solving the quasi-steady compressible Rayleigh equation 

p" - ^p' - [a 2 (l - M 2 0 ) + /? 2 ]p = 0, (2.15a) 


subject to 


p — y 0 [or outgoing waves] as y — » oo, 

2 

p ~ P to 2 Ay 3 / 2 as y — ► 0+, 

O 


(2.156) 

(2.15c) 


and 


P 

A 


(«„ = , k = /” Ai{q)ig) . (2.15 i) 


(1 + W) 

In (2.15a) the prime denotes a y - derivative, while (2.15d) stems from the viscous linear 
lower-tier solution in terms of the Airy function Ai. The properties of the moderate - cooling 
system (2.15) are investigated in the next section. 
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3. INSTABILITY PROPERTIES FOR MODERATE COOLING 


Solutions of the moderate - cooling problem (2.15) were obtained (a) by computation, (b) 
analytically for small r, (c) analytically for large r, connected with relatively high frequencies, 
and these are described in turn below. In practice, some of the features found in (b), (c) 
were used to guide the computations in (a), and vice-versa. 


3(a). Computational Study 

For moderate values of the Mach number we chose the mean-flow profiles U 0 and T 0 to 
have the following representative forms: 

U 0 = 1 - e-*'\ (3.1a) 

T 0 = l-e~* L, \ (3.16) 

i.e. they have the required square - root dependencies as y — ► 0+; the decay as y — ► oo is 
of relatively little consequence, we note. Then M 0 = UqMoo/T^ 2 = M 00 ( 1 — e -5 ^) 1 / 2 , and 
the system (2.15) becomes 

f - \r'' \ ^ ~ (“’O - - o'*") + ?)P = (32a) 

with the boundary conditions 


and 


p — > 0 (or outgoing waves) as y — > oo, 
2 

p ~ P rc^Ay 3 ^ 2 as y 0+, 

o 

P (ia)- 1 / 3 Ai'(£ 0 ) 


(3.26) 

(3.2c) 

(3. 2d) 


A (l + /?Va 2 ) <to) ’ 
where £ 0 and k are defined by (2.15d). For y « 1 we may obtain two independent solutions 
of (3.2a), with multiplicative constants a\, 02 , 

1 ^ n 2 M 2 

JS = a,{l + (/? + - V + “ ! (1 + + a 1 ) + + • • ■} 


+<‘Af ,a - + ijf 6 " - +-}, f°r V « I- (33) 

One numerical method used to solve the eigenvalue problem (3.2) is to start the integration 
procedure at a small value of y = 8, say. Then for non-neutral solutions, for example, we 
can fix /3 and guess a, complex, and find two independent solutions for p corresponding 
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to (ai,aj) = (1,0), (0,1) using a standard 4th order, 4 stage Runge-Kutta scheme. The 
integration is continued to a large value of y, say, where p ~exp[-(o J (l -Ml) + P)'l'y] 
and the two solutions p x and p 2 are combined in the form 

2 

p = Ppi - -tcl 2 p 2 , (3.4) 

where P is defined by (3. 2d). Without loss of generality we may choose A = 1. This system 
is iterated until the condition 

p'{y oc) + ^p(y 00 ) = 0 (3.5) 

is satisfied, where A = (a 2 (l — M^J + Z? 2 ) 1 / 2 with Re(X) > 0. In another method, to avoid the 
possibility of exponential growth in the numerical solution on the interval (5, Voo), the solution 
of (3.2a) can be found first from y — 6 to y = y 0 (y 0 = 1, for instance), choosing A = 1 and 
using the asymptotic behavior of (3.3) as initial conditions. Additionally a solution of (3.2a) 
may be found from y = y^ to y = y 0 using the exponential decaying behavior at y = y OQ) 
and choosing the constant multiplying the solution so that the two solutions are continuous 
at y = y 0 . The condition that p and p> are continuous at y = y 0 is sufficient to determine 
the complex eigenvalue a for fixed /3, r and fi. In some instances, to aid convergence, 
a further method was used where it was found useful to map (0,oo) onto (0, 1) by use of 
the transformation y = tan(7rs/2). The method of solution in this case can be either of the 
two methods described above, since the asymptotic behavior of the transformed equation of 
(3.2a) as s — ► 0+ and 3 — > oo is unchanged (simply replace y by tan(7T3/2)). Each method 
described above was found to give results consistent with the other two. Each scheme was 
used with different values of 6 and the step length for the Runge-Kutta integration scheme 
to test for convergence. 

Neutral and unstable solutions of (3.2) were found by using the numerical methods de- 
scribed above and are shown mainly in Figure 2; see also in Figures 3,4 later. The behavior 
of the solutions was followed as the (real) frequency fi was increased. The neutral solutions 
were obtained for two-dimensional and three-dimensional disturbances. It was found useful 
to plot the neutral curves for fixed values of /3/a. Figure 2(a) shows the neutral values of r 
and for M 0 0 = 0.5 and ft l a — 0, 1, 2, 4. We see that r decreases as fi increases. The neutral 
solutions for supersonic flows differ depending on whether (3/a — ( M ^ — l)s is positive or 
negative. (The situation with /3/a = (M^, — l)a is not addressed here.) Figure 2(b) shows 
neutral values of r and f2 for Afoo = 2 and /3/a = 1.75, 1.5, 1.25, 1,0. For /3/a > (M£, — 1)5" 
the solution for p is real and decays exponentially as y — > oo. For /3/a < i M lo - 1)* we 
allow outgoing traveling waves and so p is complex and has a purely oscillatory behavior as 
y — > oo. These two different types of solution produce different forms of the neutral curves 
as shown in Figure 2(b). For /3/a > (M^ — 1)? (/3/a = 1.75) there is no solution for large 
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values of both r and ft. For 0/a < (M^ — 1)? there exists a cut-off value of ft for each 
value of 0/a below which there are no solutions. As ft increases from the cut-off value, 
on the upper curve, r reaches a local maximum and then decreases to a local minimum 
as ft increases further. Beyond the local minimum the curve has no more turning points 
and r continues to increase as ft increases. We note that the upper curve for the 2D case 
0/a — 0 does have this behavior but because of the scale of Figure 2(b) it is not immediately 
apparent. The lower curves become almost coincident as ft increases, including the cases 
/3/a > ( M — 1)? for which r continues to grow as ft becomes increasingly small. Thus we 
see that there exist neutral (and hence unstable) solutions of (3.2), over a range of r, for all 
values of ft for (3 /a > (M£, — l)s or for any subsonic case, but there is a cut-off value of ft 
for each value of /3/a for which /3/a < (M^, — 1)? where no neutral (or unstable) solutions 
exist for ft below this value. Figure 2(c) shows neutral values of a and ft for M = 2 and 
/3/a = 1.5, 1.25,1,0. The curves for (3/a > (M£, — 1)? for the supersonic case and for any 
subsonic case are also shown. They are only present on the upper curve where the different 
solutions are indistinguishable. For /3/a > (M^, — 1)? and also for M < 1 the pressure 
p is real, hence, from (3. 2d) P must be real also. This implies the well known result that 
£ 0 = —2.298 i? and i~* Ai'(£o)/K,(£o) — 1.001. The definition of £o from (2.15d) determines 
a for these cases. We find that a = (ft/2.298)*. Hence, for the cases where this solution is 
valid, a is independent of Moo, 0 and r. The two branches of the neutral curves shown in 
Figure 2(c) correspond to the asymptotic solutions described later in sections 3(b) and 3(c). 

The effects of increasing the cooling of the surface were investigated by decreasing r for 
two-dimensional, three-dimensional, subsonic and supersonic (for moderate values of Moo) 
disturbances. The trends were found to be the same in all of these cases, namely that 
decreasing r increases the growth rate of the unstable solution. This is shown, for example, 
in Figure 2(d) for Moo = 0-5,0 = 1 and r = 0.1,0.05,0.02, where the growth rate —Im(a) 
is plotted against ft. We observe that the growth rate increases significantly as r decreases 
(i.e. when the wall is cooled more). This suggests a large destabilization of the basic flow 
as the wall cooling is increased. The results for Moo = 0.5 and /3 = 0 for r = 0.1, 0.05,0.02 
are graphically identical to the respective results for 0 = 1 shown in Figure 2(d). The real 
part of a is very large for these solutions for small values of t. Figure 2(e) shows the growth 
rate — 7m(a) as a function of ft for M 0 0 = 2, 0 = 1 and r = 0.1,0.05,0.02. The solutions 
have the same behavior as those shown in Figure 2(d) for M^, = 0.5 with the growth rates 
appearing to be slightly smaller for M^ = 2 compared to those for M^ = 0.5. Again, the 
corresponding solutions for = 2 and 0 = 0 are virtually the same as those for M x = 2 
and 0 = 1 for the values of r shown in Figure 2(e). 

For large values of r the solutions approached those of Smith (1987) as expected j»ee 
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also (c) below]. For subsonic flows as r increases the solutions for Im(a ) have the same 
behavior as those shown in Figure 2(d) with the growth rates continuing to decrease as r 
increases. This behavior is found to hold for very large values of r. This is illustrated in 
Figure 2(f) where — Im(a) is plotted against for = 0.5, (3 = 0 and r = 0.4, 0.5,1. 

If the corresponding curves for /? = 1 were to be added to Figure 2(f) they would be 
indistinguishable. Figure 2(g) shows —Im(a) as a function of for M & = 0.5 and (3 = 0 
and (3=1 for r = 10, 20. The solutions for (3 = 1 and (3 = 0 are now distinct, with the 
solutions for (3 = 0 having larger growth rates, the difference increasing as r increases. We 
note that as 0 increases from zero the three-dimensional mode occurs first. 

For supersonic flows, as expected from the neutral curves shown in Figure 2(b), the 
structure of the eigenvalue a. changes as r increases. An example of this is shown in Figure 
2(h) where — 7m(a) is plotted against f l for M «, = 2 and /? = 0 for r = 0.4, 0.5, 1. For 
0 < r < 0.4 there is one neutral point at a small value of fi. For 0.4 < r < 1 there are 3 
neutral points and for r > 1 only the neutral point at the largest value of f2 persists. Figure 
2(i) shows —Im(a) as a function of for larger values of r, namely for = 2, /? = 0 and 
t = 5, 10, 17. For the range of v shown in Figures 2(h) and 2(i) the solutions for (3 = 1 are 
very close to those shown for /? = 0. As r increases the growth rate continues to decline. 

To summarize the results of the computational study for moderate values of we 
can make the following general statements. As the amount of wall cooling is increased 
(i.e. as r is decreased) the growth rates increase for two-dimensional and three-dimensional 
perturbations, irrespective of whether the flow is subsonic or supersonic. Moreover, as r 
is increased, for subsonic flows no significant three-dimensional effects are noticeable until 
r > 10, with the growth rates here being larger for two-dimensional perturbations. On the 
other hand, for supersonic flows with (3/a < ( M — 1)5 there is very little difference between 
the solutions for two-dimensional and three-dimensional flows. Again, for a fixed value of r 
the growth rates for subsonic flow are larger than those for supersonic flow. 

We also note that, here and below, multiple modes may be present but these have not 
been detected as yet. 

For large values of M ^ the mean flow becomes different in structure. In particular, from 
the Howarth-Dorodnitsyn transformation, To ~ so that y ~ M^,, while Uq ~ 0(1). We 
do not present solutions for large Mach number but below is an outline of a representative 
problem to be solved in this case. For large we write 

(3.6) 

and then 

[00, To, Mo, a] = [0o(F), MlT 0 (Y), M 0 (F), M^ 2 a\, (3.7) 
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instead of (3.1a,b), with the over-barred variables now being typically of order one. The 
equation (2.15) therefore becomes 

p r7 -?^p r -[a ! (l-M;) + 3 I ]p = 0, (3.8) 

for 0 < F < F TO , where the scaled boundary-layer edge F^, is finite in this regime, M 0 = 

and /3 = M^ 2 p. The boundary conditions near the cooled surface likewise remain 

intact, to leading order, with z, P, A also scaled appropriately with respect to M^, taking 

account also of the extra factors such as Ci, c 2 involved in the scales of section 2. In particular, 

6 

the typical wall-cooling factor r required now increases like MJL- A logarithmically thin extra 
layer near Y = Y adjusts the solution from T 0 being 0(M£J (as above) to being unity as 
required in the free stream, and so the profiles Uo, T 0 are to satisfy the conditions Uq — * 1, 
T 0 -> 0 as Y -> F«,- and U 0 <x F*, T 0 oc F^ as F 0+. 

(b) Small r: towards severe cooling 

In line with the computational results, it is suggested that for small values of r the 
disturbance y - scale contracts. So U 0 , T 0 are approximately like y 1 ^ 2 , M 0 is small of order 
jjf 1 / 4 for Moo of 0(1), cf. later, and therefore (2.15a) suggests that a is 0(|T 1 ), as is ft. 
The unknown y-scale is then fixed in terms of r by comparing the typical size of P/A from 
(2.15c) with that inferred from (2.15d) , given that is expected to remain typically 0(1) 
for instability. The implied small-r response therefore has 

p = p(y ) + •••,. y = T 6 / 5 y> ( 3 - 9a ) 

(a, p) = r- 6 / 5 (d, P) + • • • , ft = r-“/ 6 ft + • ■ ■ , (3.96) 

p = /> + •••, A — A -+-•••. (3.9c) 

This reduces (2.15) to the equation 

m-f?-(“" + /5 1 )p=0. (3.10a) 

21/ 

for p, subject to 

p( 00 ) = 0, p ~ P — -a 2 Ay 3 ? 2 as y — » 0-|-, (3.106, c) 

O 

and 

P/A = (1 + P'/a'r'iiay^At'tfoyKtio), (3.10*0 

where £ 0 = -i 1/3 ft/d 2/3 . Here the solution of (3.10a,b) can be obtained analytically, 

P = diy 3/4 A|(y) with y = (a 2 + P 2 ) 1/2 y, (3.11a) 
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di is a constant, K 3/4 is the Bessel function in the standard notation, and then (3.10c) 

requires that 

2ttq. 2 

(r(i))2(«2+ £2)3/4 ’ 

from the behavior of the Bessel function at small y. So the coupling with (3.10d) leads to 
the reduced eigenrelation 



P/A 


^(£ 0 ) _ 2?r 

*«<>) ~~ (r(j )) 2 


(«) I/3 ( 


*2 

a 




(3.11c) 


between a, Cl. The solution for i2e(a) and Im(a ) from (3.11c) are plotted in Figures 3(a) 
and 3(b), respectively, for (3 = 1 and r = 0.1,0.05,0.02. These are seen to agree well with 
the full computations at small values of t; see also in Figures 2(d) and 2(e). In order to 
illustrate this agreement Figures 3(c) and 3(d) show Re(a ) and Im(a ) for r = 0.05, (3 = 0, 
for Moo = 0.5 and Moo — 2 from the full computations along with the result of (3.11c). The 
agreement is seen to be very close when Cl is not very large. 

It is interesting that the small-T response here is independent of the Mach number, as 
might be expected for disturbances focussed near the surface. This is because the effective 
Mach number M 0 is small, so that the contribution M 2 in (2.15a) is negligible, although the 
instabilities that are found to exist do remain temperature- and velocity-driven through the 
ratio M 0 y/M 0 operating in (2.15a). The property that these extra-cooled- surface instabilities 
exist at any Mach number, and for both 2D and 3D modes, contrasts with the situation for 
the less cooled surfaces corresponding to large r or to the earlier regime of (2.1), where, in 
the supersonic range, instabilities occur only for sufficiently oblique modes directed outside 
the wave-Mach-cone: see (c) below. Thus the cooling now destabilizes a range of previously 
stable viscous modes. Moreover, the de-stabilization is very marked, for all the viscous 
modes, in the sense that the spatial growth rates produced are much larger than encountered 
before. Indeed, the growth rates appear to exceed those of purely inviscid modes (see also 
section 5). The small-T response leads on to the new structure investigated in section 4. 


(c) Large r and the cut-off 

The clue to the behavior at large r, associated with less cooled surfaces, is found in 
high-frequency properties. Thus when Cl is large, at 0(1) values of r first, the outer solution 
of (2.15a) takes the form 

[p, P, A, a, 1 3 , Mo] = [p, P, IT 1 / 3 i, Cl 2 ! 3 a, Cl 2 ' 3 /. 3, 0(CI~^ y 1/4 )] + •■•, y = n~ 2 ' 3 y, (3.12a) 
concentrated again near the wall. So p satisfies the equivalent of (3.10a), subject again to 
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(3.10b,c). Here however the sublayer relation is simply 

P/A = a -1 (l + fi/a 2 )- 1 , (3.126) 

as |£ 0 | is large for these high frequencies. Hence we have the result 

° = ( r (j)) 4/3 /(2iTT) 2/3 , (3.12c) 

for the 2D example [the 3D version has an extra factor (1 -f /? 2 /a 2 ) -1 / 6 on the right-hand 
side], giving neutral stability at leading order. This result is essentially an inviscid one since 
the asymptote — £o for Ai /k in (3.10d), used in (3.12b), corresponds to an inviscid limit, 
and, like that in (b) above, the present result applies for all Mach numbers. 

Closer inspection then shows that the high-frequency structure and results (3.12b,c) 
continue to hold for all r, provided that Q >> t - 4 / 5 and D >> r, the former to keep |£ 0 | 
large and the latter to keep a large (and y small), as supposed above. The former restriction 
is associated more with decreasing r since (b) comes into play exactly when 0 reaches the 
order r -4/5 and thus provides a match. The other restriction is connected with the large-r 
behavior below. Also, we note that higher-order contributions in inverse powers of D bring 
in the second term in Ai' / k, which is a viscous and de-stabilizing effect, such that there 
exists instability for any value of r, at sufficiently high frequencies. 

To complete the picture at large t, then, we move on to consider Q values of order r 
(Figure 4). The properties are mainly inviscid again since |£ 0 | is large, in fact comparable 
with 0, as 

r ~ fir with a,/3 = 0(1), y = 0(1), (3.13a) 

to leading order. The dominant equation to solve is therefore exactly (2.15a), with the 
constraint (2.15b) still, but (2.15c) becomes 

p ~ P[1 - ^fa(a 2 + 0 2 )^ /2 ) as y — > 0+, (3.136) 

because the sublayer response gives here DA = (a 2 + j3 2 )Pfa. The computed solution of 
(2.15a,b) with (3.13b) is presented in Figure 4(a) along with comparisons with the compu- 
tations of section 3(a) and the solution from (3.12). At small r, corresponding to reduced r 
and/or larger D values, there is a match with (3.12) as expected, for any Mach number. Fig- 
ures 4(b,c,d) show solutions from the computations of section 3(a) for large r with M ^ — 2 
and /3 = 0. In Figures 4(e,f,g) we show f as a function of fi along the neutral curves for 
Mqo = 2 and /? = 0,1, 1.5, as shown in Figure 2(b). We see that for the branch r >> 1 
and D >> 1 f decreases as D increases so there will be a cut-off value of f above which no 
unstable solutions exist. As f increases the solution becomes Mach-number dependent, and 
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this is especially so at large f where the subsonic realm remains always unstable but the 
supersonic realm becomes stable within the wave- Mach-cone /3 < a(M£, — l) 1 / 2 , including 
the 2D case. The cut-off here, for the less oblique supersonic modes, is profile-dependent 
but occurs at an 0(1) crossover or neutral value of f. That yields the predictions (in line 
with the computations just mentioned and with those of Figure 2) 


Cl — neutral oc r 
Cl — max oc r 


as t — ► oo, 


(3.14a, b ) 


for Moo > 1 and waves directed inside the wave-Mach-cone. This behavior is illustrated in 
Figure 4(h) where Cl — max is plotted as a function of r for f3 = 0 and M ^ =0.5. The 
frequency = fi — max is that at which the maximum growth rate occurs. In contrast, 
analysis for the more oblique supersonic modes and for the subsonic ones shows that at large 
r the outer solution develops a double structure, leading to the typical response 


a ~ 


(1 - Mlf'Vr 1 ' 2 , 


(3.15) 


for the 2D subsonic case for instance. This confirms the suggestion of the computations, 
that the result (3.14b) in fact holds for all the supersonic and subsonic modes, whereas the 
subsonic and more oblique supersonic ones retain instability at large f (i.e. Cl << r), as 
opposed to (3.14a). 

The double structure just mentioned leads to the final phase for the latter modes, namely 
when Cl falls to 0 (t~ 1/f2 ). Then |£ 0 | decreases to 0(1) since 

(a, /3) ~ r _3/4 (a*,/9,) with ~ r~ l ^ 2 Cl,, (3.16) 


from (3.15), with a,,/3,, Cl, of order unity, so the fully viscous sublayer relation (2.15d) is re- 
installed. Moreover, the outer double structure produces the two y-scales 0(r 3//4 ) and 0(1), 
the former giving a potential-flow character and the latter a profile-dependent displacement 
effect, leading to the relation P = a,Ar~ 3 ^* (1 for subsonic 2D modes for instance. 

Coupling this with (2.15d), then, we obtain the eigenrelation 


T&ki™-) 1 '* = (1 

At (< 0 ) a. 


(3.17) 


between a,, Cl,, again for the 2D subsonic case as an example. There are two main points 
here. First, the relation (3.17), which agrees with (3.15) at large Cl, as required, produces a 
maxi m um growth rate at an 0(1) value (approx. 4) of Cl, and neutral stability at a smaller 
0(1) value (Cl, approx. 2.3), thus giving rise to the predictions 

Cl — neutral ~ r -1 / 2 


as r 


oo. 


fl — max 


r-1/2 


(3.18a, b ) 
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These again appear to be consistent with the computations at large r, as shown in Figure 4(i) 
where fl-max is shown as a function of r for /? = 0 and M , = 2. Second, the entire structure 
in this regime, and the results (3.16) - (3.18), match directly with the earlier triple-deck case 
(see (2.1)). For instance, the wavelength 0(Re~ l ^ 2 a~ 1 ) increases to <9(Re -3 / 8 ) when T w is 
raised to 0(1), since then T w oc r 1 / 2 -► 0(i2e 1/12 ) and a — > 0(Re -1 / 8 ) from (3.16). Similar 
matching holds with regard to (3.17), (3.18), for 2D and 3D subsonic or supersonic modes 
(Smith 1987). 

The splitting of the neutral curves for the different cases, and the production of two local 
maxima (and mounds of instability) in certain cases, via (3.14), (3.18), reflects the emergence 
of two major types of instability when the surface temperature increases. These are the 
viscous-inviscid triple-deck type (of (2.1) ff), with which (3.18) connects [see comments 
above], and the inviscid compressible- Rayleigh type, which are implied by (3.14). In the 
latter case, raising T w back to 0(1) corresponds to T w increasing to 0{Re 1 ^ 12 ), hence r 
to 0(Re 1 / 6 ), and so the fi’s in (3.14) become 0(i2e 1 ^ 6 ) also, meaning that the time scale 
t becomes 0(Re~ 1 / 2 ), from (2.8). This t-scale, the x-,y-scales of order Re -1 / 2 , and the 
increased wavespeed of order unity, all coincide with the compressible Rayleigh scales : see 
also section 5. 

Finally here, a combination of the small-T features in (b) above and those for large t in 
(3.14) indicates the existence of an 0(1) value of the effective surface temperature at which 
D-neutral is a minimum, for the supersonic 2D and less inclined 3D waves. This is as shown 
in the figures. The corresponding fi value in the 2D case for instance represents a critical 
value of the scaled frequency, below which 2D modes remain stable, in supersonic flow, and 
above which they are unstable for a range of surface temperatures. 

4. THE SEVERE-COOLING REGIME 

This regime has much reduced surface temperatures and its structure is inferred from the 
properties found in section 3(b). Thus for small r the characteristic thickness of the outer 
tier deceases like r 6 / 5 , as does the streamwise length scale, while the buffer thickness which 
is proportional to T 2 decreases like r in view of (2.12a). Along with this, the representative 
viscous-tier thickness falls as a~ l / 3 T w because of the scaling in (2.7) and the definition of 
The unsealed thicknesses are therefore of the orders Re _1 / 2 T 6 / 5 , Re~ 2 / 3 r, Re~ 3 / 4 T 9 / 10 
respectively. 

The severe-cooling regime arises when any two of the last three thicknesses mentioned be- 
come coincident. For instance if the buffer thickness coincides with the viscous-tier thickness 
then r Re 5 / 6 , giving 

T w ~ Re- 1 / 2 (4.1) 
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as the typical surface temperature here. At that stage the other thickness above also becomes 
comparable with the other two, and so a single-tier structure is expected to apply for severe 
surface cooling. In effect the original main tier and buffer layer now enter the viscous 
sublayer. The corresponding streamwise length scale is 0(Re~ 2 / 2 ) which is relatively tiny 
and coincides with the inner-tier thickness where y ~ Re~ 3 ^ 2 and viscous compressible effects 
are felt. The corresponding velocity, pressure, temperature, density and viscosity scales can 
be worked out likewise. The resultant governing equations for the severely-cooled flow are 
in effect the unsteady compressible Navier-Stokes equations, 

Pt + {pu) x + (pv)y + (pw) z = 0, (4.2a) 

p[d t + ud x 4- vd y + wd z }u = -p x + d x [2fj.u x + (// - ^/x)V.VJ + d y [p(u v + u x )] + d z [p(w x + u z )] } 

(4.26) 

p[d t +ud x +vd v +wd z ]v - -Py+dy[2p,v y +(p,'--fj,)'V.V]+d z [n(v z +w y )]+d x [p.(u v +v x )], (4.2 c) 

p[d t + ud x + vdy + wd t ]w = -p z + d z [2pw z + (/z' - ^p)V.Y] + d x [p(w x + u,)) + d v [p(v z + w v )], 

( 4.2 d) 

p[H t + uH x + vH y + wH z ] = V.(k gradT + ru), (4.2e) 

a = + V 1 + t» s ) + V - fi = CT, p = KpT, (4.2/, g, h) 

L p 

where H is the enthalpy, r is the stress tensor and a is the Prandtl number. Here again 
the nonlinear interaction is set up first, with small perturbations then being governed by 
the linearized version. The basic undisturbed flow has u = T^cacJ -1 [( 1 + 2c 1 y) 1 ^ 2 — 1] and 
T = 7^(1 + 2c 1 j/) 1 / 2 , as in (2.3) in effect, and the boundary conditions on (4.2) include 

U ~ V = w = 0, T — T W} at y = 0, (4.3a) 

u ~ T w —[(l + 2 c iy ) 1/2 - 1], T ~ T w ( 1 + 2c 1 y) 1 ^ 2 , p — ► 0, as y — > oo, (4.36) 

Ci 

where T w is the scaled surface temperature. 

The main difference from the previous cooling regime is that the viscous Navier-Stokes 
motion is now fully compressible (and single-structured), contrasting with the quasi-incom- 
pressible sublayer properties holding before. The scales involved are extreme but they can 
be verified from an order-of-magnitude argument. Further, the present flow structure applies 
for any Mach number, in line with the closeness of the Navier-Stokes structure to the surface; 
and both the spatial and temporal growth rates induced for small perturbations are an order- 
of-magnitude larger than those found previously for inviscid modes, because of the tiny scales 
involved here. 
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This new predominance of the viscous modes is due physically to the high surface shear 
stresses and surface heat transfer rates that the cooling of the surface provokes. The 
same comment applies to section 3(b). The precise behavior is of course dependent on 
the temperature-viscosity law in operation (unlike the law-independent triple-deck case of 
(2.2)), but the same qualitative behavior seems likely to hold in general, for the same physical 
reasons. Linearized or nonlinear solutions in this regime could therefore be of some interest, 
as T w is varied. 

5. ON INVISCID MODES ON COOLED SURFACES 

The pure inviscid stability of compressible boundary layers is controlled by the compress- 
ible Rayleigh equation 

p"-=?-[a’(l-]B*)+0 l ]ii = O, (5.1a) 

for the pressure perturbation, with the length scales as for the main tier in section 2 but 
now the waves peed c is 0(1), so that 

M = (U 0 - c^To 1 ' 3 , (5.1 b) 

cf. (2.15a). The boundary conditions are to match with the undisturbed incident free stream 
and to satisfy the tangential-flow surface condition, 

p — ► 0 [or outgoing waves] as y — > oo, (5.1c) 

p* = 0 at y = 0, (5. Id) 

respectively. Computational solutions of (5.1) are given by Mack, Malik and others (see 
references in §1) including some investigations for cooled surfaces. The investigations trace 
the effects of gradual cooling on the originally un-cooled inviscid modes and find some inter- 
esting features. There is no significant change of scale however. The present investigation 
focuses mostly on new inviscid modes involving substantial scale changes and, in particular, 
much increased growth rates. 

The main factor for cooled surfaces where T w is small is the development of the 0(T 2 ) 
sublayer described by (2.3), associated with high wall shear and temperature gradient. The 
previous investigations mentioned above neglect the sublayer in a sense and seek solutions 
of (5.1) on the scale y of 0(1). Here, in contrast, we find significant effects arising when the 
sublayer plays an active role. We concentrate now on the two-dimensional case. 

The first type addressed has c » TjJ 2 and a 2 c 2 ~ T~ 3 . This is suggested by the form 
of M in the sublayer, namely 0(T 'l/ 2 ) — c times 0(T~ 1 ^ 2 ) from (5.1b) with (2.3a), and by 
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(5.2a) 


the balance of terms in (5.1a). So M is large in the sublayer and (5.1a) reduces to 

k 2 p 


—II I P , - r _ n 
p + rr + tt/2 ~ °» 


2y 


y 


where y = y + (2ci) 1 , the parameter k = acM 00 (2c 1 )~ 1 ^ 4 T^/ 2 is 0(1) and c x is the heat- 
transfer factor defined in (2.3b). Here (5.2a) becomes Airy’s equation in terms of the variable 

A 1/2 

s — y and so the constraints (5.1c,d) lead to the eigenrelation 


-2iAi(-(4fc 2 ) 1/3 ) + Hi(-(4fc 2 ) 1/3 ) = 0 


(5.26) 


for k. This has no real roots. All the (complex) roots give damped waves in fact, but this 
type of mode serves to suggest the new scales for the second type of surface-cooled mode. 
Here c reduces to 0(T^/ 2 ), and a ~ T~ 2 , as implied by the first type, and this brings in a 
new effect since now ~M is of order unity. Thus if c = T* /2 c and a = T~ 2 a then JT oc y _1/ * 
to leading order and (5.1a) becomes 


p" 


, f , -2r 

T* — A T OL A 

2 y (2 Cl yy/ 2 


- l]p = 0. 


(5.3a) 


A 1/2 

For a real, the substitution y = c 2 M^,/(8c 1 ) 1 / 2 + 5/ (2a 1 / 2 ) converts (5.3a) to the parabolic 
cylinder equation 

s 2 - v 

Pr,-(j-P)p = 0, (5.36) 

where ft = ac 4 M^,/(8ci) and the boundary conditions are p(oo) = 0 ,dp/ds = 0 at J = 
(2 — c 2 M% 0 )a 1 / 2 /(2c 1 ) 1 / 2 . The solution is therefore the parabolic cylinder function p = U(a, 5) 
with a = —(3 (see Abramowitz & Stegun (1964)), for the outer constraint, while the inner 
constraint requires that s = x n where x n are the roots of 


U'(a,x n ) = 0. (5.3c) 

The above reference shows that there do exist some roots with a being negative, (3 positive, 
as required, occurring for positive x n . For those roots the values of a, c then follow from the 
relations 

(2 - c^Ja^/pc,) 1 ' 1 = *„, (5,4a) 

ac'M'J{2c x ) = 0. (5.46) 

The resulting wave speeds c are all real, given by 

3 1/J = 03>/2 + i n /2)(2c 1 )'/ 1 , M„c = 2(2 + xJp'l*)-'/\ (5.5a, 6) 
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where the first equation fixes a, the second c. Hence this new type, which is a temperature- 
driven mode because of the form of M, produces new inviscid neutral modes. 

The third case has c decreased still further, to 0(T W ), which reinstates the influence of 
the velocity profile Uq as well as the temperature at leading order. For, if c = T w c, a — T~ 2 a, 
then M = T£/ 3 M is small and the controlling equation becomes 

>=. 9 

p"~ - ai 2 p = 0, (5.6a) 


with 

M = [{y^ 2 - 1)2A X - c]Mooyi~ 1/4 . (5.65) 

(Here y x = 2 c x y, a = 2ciai, 2c x \ x = c 2 ). This can be recast, by setting y x 1 ^ 2 = /2(1 + £), fi = 
1 + cj 2Ai, in the form 

(Pti - % - + !) 2 P = 0. (5.6c) 

with 

p(oo) = 0, pf = 0 at ( = /i -1 — 1, (5.6d) 

where (f> = 2a x fi 2 acts as the eigenvalue. This third type can be shown to match with the 

second type, at large c, ft, and similar matching occurs between the first and second types. 

Hence again in this regime neutral modes are produced. 

The fourth type of cooled-surface inviscid mode to consider has a structure suggested by 

the work of section 3(c). Here c is reduced to 0(T 2 ) but a is 0(1), and the solution of (5.1) 

takes on a three-tier form. In the bulk of the boundary layer where y is 0(1), therefore, 

(5.1a) applies but with c replaced by zero in the definition of M in (5.1b), i.e. (2.15a) holds 

1 /2 

again with M 0 = UoMoo/T 0 , and now 

p ~ P — ^a 2 A**!/ 3 / 2 as y — ♦ 0+, (5-7) 

O 


where the relation between the constants A** and P has to be found. The next tier, closer 
to the surface, has y = T 2 y along with (2.3), so that to leading order 


M = T'J 2 M{y ) : 


M= c 2[(l + 2ci p) 1/2 - 1 ]Mqo 

Ci(H-2c 1 y) 1 / 4 


(5.8a) 


Now 

P = P + T 2 Pl + • • ■ 

and consequently the first two terms in (5.1a) dominate, giving on integration 

dpi aM«* C ; /a J a 

dy ^M^c 2 ’ 


(5,85) 


(5.6c) 
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where the coefficient on the right side is chosen to match with (5.7) and we normalize such 
that pi(0) = 0. Hence dpi/dy decreases like y 3 at small y, which leads on to the lowest tier 
even closer to the surface where f = T 3 y with y of 0(1). There 

f=P + T*p 1 + ---, (5.9a) 


and the small wavespeed c = T 2 c takes effect since 

M = T^] 2 M + • • ■ with M = (c 2 y — c)M OOJ 


(5.96) 


from (5.1b) with (2.3a). Hence the governing equation, from (5.1a), is 

= a 2 P. 


d 2 pi 2 dM dpi 2 


dy 2 M dy dy 
The solution for dpi/dy merging with (5.8b,c) at large y is then 


dpi 

dy 


a 2 P , . a 2 A**c 3 i /2 f _ _ v 

-{c 2 y - c) - 21/2c — {c 2 y - c ) 




C 2 ^ ’ ^2 

and so the required surface condition (5. Id) is satisfied at y = 0 if 

2 1 / 2 c 2 

A = t4 r 

c 


(5.9c) 


(5.9d) 


(5.10) 


T his provides the relation for A**, to be substituted into (5.7). The resulting problem for p, 
that is, (2.15a, b) with (5.7), (5.10), is equivalent to the one derived and solved numerically 
in section 3(c): see (3.13b), where the parameter fa 3 involved for the two-dimensional case 
is in agreement with that inferred from (5.7), (5.10), with account taken of the different 
temporal scaling in (2.8). Hence the stability properties in subsonic and supersonic flow, for 
this fourth type, are exactly as given in Figure 4 with the appropriate change of notation. 
Like all the types described in this section, this new type appears to be quite distinct from 
any cooled modes found previously. 


6. FURTHER COMMENTS 

The main conclusion from the present theoretical research [in particular, §§3,5] is that 
moderate surface cooling can enhance the TS viscous-mode activity to the extent that the 
spatial growth of the viscous modes can become comparable with, or even exceed, that of 
inviscid ones, thus reversing the normal order. In addition, the cooling can completely desta- 
bilize otherwise stable viscous modes and inviscid modes, at any Mach number, including the 
production of radiating waves in the supersonic realm. In those two senses, then, cooling of 
the surface seems very “dangerous” with regard to transition in the compressible boundary 
layer. 
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The TS growth, when it exceeds that of the pure inviscid modes studied previously, 
occurs for disturbances which are confined relatively close to the surface, as opposed to the 
full inviscid ones stretching all the way across the boundary layer. In particular, for all 
the new modes, the high-heat-transfer sublayer plays a substantial active part in enhancing 
the viscous growth, through the correspondingly high shear-stress production at the surface. 
It could be for that reason, i.e. the appearance of multi-structured responses within the 
boundary layer, that previous theoretical, computational and experimental studies of surface- 
cooling effects appear to have not detected the new mode properties, although there are other 
explanations which suggest themselves such as considerable finite-Reynolds-number effects; 
see also below. Experiments and computations focusing on the possibility of such strong near- 
surface behavior remain to be done, it seems, and they could be of much interest scientifically 
and for practical applications. There is some measure of agreement, nevertheless, between 
the present theory and the Lysenko and Maslov (1984) experiments. Thus the theory predicts 
that as the temperature ratio T w falls from the non-cooled regime to the moderate-cooling 
regime, for first modes at fixed Reynolds number, the typical streamwise wavenumber, a D 
say, and the typical frequency, u D say, increase as T~ 3/2 , T~ x respectively: see §2. So for the 
1 to 0.76 decrease in T w examined in Lysenko and Maslov’s Figure 2 the predicted increase in 
u D is by a factor 1.316, which is not inconsistent with the measured increase in the neutral 
value of w D shown in that figure. (See also comparisons for the insulated case in Smith 
1987). Similar trends for decreasing T w are observed in their Figure 5. Again, Lysenko and 
Maslov’s Figure 4 has a case of a 1 to 0.79 decrease in T w , for which the present theory 
implies an increase by a factor (0.79) 3 / 2 , i.e, 1.398, in the corresponding spatial growth 
rates — af\ This also is not inconsistent with the measured growth curves in their Figure 4. 
Hence as far as their Figure 4 is concerned there appears to be some agreement both on the 
delay of instability and on its increased growth rates, with cooling of the surface, although 
their interpretation (see their figures) is that the growth rates in fact decrease. The latter 
would seem to be due to extra normalization factors used, whereas the theory indicates that 
the effects of the temperature ratio T w should be the overwhelming ones eventually, as T w 
falls. 

A further note of caution should be added here, concerning nonparallel-flow effects. We 
recall that, strictly, supersonic first modes within the wave-Mach cone are stable and gov- 
erned at leading order by nonparallel-flow effects, for T w of order unity (Smith 1987). There- 
fore the experimental results (of Lysenko and Maslov (1984), and others) inside the wave- 
Mach cone could well be influenced to a significant extent by nonparallelism rather than being 
true instability measurements, and they may in turn miss the new supersonic instabilities 
predicted for cooled surfaces. 
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The corresponding nonlinear processes at higher disturbance amplitudes which are noted 
earlier in the text should also be of much interest with respect to transition in the cooled 
boundary layer. For example, the nonlinear break-ups described by Smith (1988) then be- 
come attainable within finite time; and forms of vortex-wave interaction can also be induced, 
including those of Smith & Walton (1989) and Hall & Smith (1989). As well as the nonlinear 
problems mentioned in the text, that inferred from section 3(c), for reduced surface tem- 
peratures, is governed by the three-dimensional unsteady boundary-layer equations coupled 
with the pressure-displacement law which is given in transformed terms by (3.5b). This can 
be expressed as an integral relation between the unknown pressure and displacement, in two- 
dimensions or three-dimensions. Further, the finite time for nonlinear break-up referred to 
above becomes shorter as the surface is cooled, suggesting that intermittency (see references 
above) may arise that much sooner on a cooled surface. 

There are many other regimes of concern. These include the hypersonic realm with its 
strong vorticity mode (Brown & Smith 1989); the transonic realm and its peculiar instability 
features (Bowles Sz Smith 1989); cross-flow effects; upper-branch properties, e.g. between 
the mounds of instability discussed in section 3(c); other temperature- viscosity laws, given 
that the Chapman assumption is vital in the new near-surface layers, although again other 
laws of concern seem likely to induce broadly similar behavior; and the upstream-influence 
analogue which is currently being considered. These additional aspects seem less significant 
in the long run, however, than the issues associated with nonlinearity, given the increased- 
instability properties found here for cooled-surface boundary layers. 
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2/3 b uff er tier and a 3/4 viscous lower tier. Further cooling shrinks the instability structure 
still further, as described in sections 3(b) and 4. 



Figure 2(a). 
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Figure 2(b). 



Figure 2(b). Neutral value of r versus frequency fi for M , «, = 2 and /?/a = 0, 1, 1.25, 1.5, 1.75. 
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Figure 2(c). Neutral value of a versus frequency ft for Moo = 2 and /3/a = 0, 1, 1.25, 1.5. 







Figure 2(e). 





Figure 2(f). 
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Figure 4(a). 
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Asymptotic result for x ~0(1), 

a » i , p = o 



Asymptotic result for t » 1 ,Q » 1 

p = 0, M» = 2 
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Figure 4(f). 



Figure 4(f). Neutral value of f versus frequency fi for M = 2 and (3 /a = 1. 
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Figure 4(g). 



Q. 

Figure 4(g). Neutral value of f versus frequency Cl for = 2 and /3/a = 1.5. 
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